function res=validationgraph()
S2O3s=[0.022207043 ;
0.044414085;
0.055517607;
0.066621128;
0.072172889];
Sts=0;
for i=1:length(S2O3s)
    Sts(i,1)=run_iodine_reaction_real(S2O3s(i));
    Sts(i,2)=run_iodine_reaction_real(S2O3s(i)*1.1);
    Sts(i,3)=run_iodine_reaction_real(S2O3s(i)*.9);
end
figure(1);
plot(S2O3s,Sts, '.')
res=Sts
end